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CHANDRA HIGH-RESOLUTION SPECTRUM OF THE ANOMALOUS X-RAY PULSAR 4U 0142+61 
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ABSTRACT 

We report on a 25 ks observation of the 8.7 s anomalous X-ray pulsar 4U 0142+61 with the High Energy 
Transmission Grating Spectrometer (HETGS) on the Chandra X-ray Observatory. The continuum spectrum is 
consistent with previous measurements and is well fit by an absorbed power-law + blackbody with photon index 
r = 3.3 ±0.4 and kT = 0.418 ±0.013 keV. No evidence was found for emission or absorption lines, with an upper 
limit of «50 eV on the equivalent width of broad features in the 2.5-13 A (0.95-5.0 keV) range and an upper 
limit of «10 eV on the equivalent width of narrow features in the 4.1-17.7 A (0.7-3.0 keV) range. If the source 
is a magnetar, then the absence of a proton cyclotron line strongly constrains magnetar atmosphere models and 
hence the magnetic field strength of the neutron star. We find no strong features that are indicative of cyclotron 
absorption for magnetic field strengths of (1.9-9.8)x 10 14 G. This is still consistent with the dipole field strength 
of B = 1 .3 x 10 14 G (at the polar cap) estimated from the pulsar's spindown. 

Subject headings: pulsars: individual (4U 0142+061) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

The anomalous X-ray pulsars (AXPs) have a number of 
properties that distinguish them as a class from other pulsars 
(Mereghetti & Stella 1995; van Paradijs, Taam, & van den 
Heuvel 1995). They have spin periods falling in a narrow range 
(6-12 s), luminosities of order 10 34 -10 35 erg s" 1 , and a soft X- 
ray spectrum well described by a blackbody plus a power-law 
(see Gavriil & Kaspi 2002 and references therein). They also 
undergo relatively steady spin-down, have faint or unidentified 
optical counterparts, and show no evidence of binary motion. 
Several of the AXPs are associated with supernova remnants, 
indicating that they are relatively young neutron stars. Mod- 
els to explain these properties fall into two general categories: 
accretion models and magnetar models. 

In the accretion models, the accreting matter does not neces- 
sarily originate from a binary companion. X-ray timing and op- 
tical counterpart limits allow for only very low-mass compan- 
ions, although it is difficult to reconcile such binaries with the 
supernova remnant associations because a supernova explosion 
would likely have unbound the binary. Another accretion model 
suggests that AXPs are isolated neutron stars accreting from a 
fallback disk of material from the supernova explosion or from 
the remains of a Thorne-Zytkow object (e.g., van Paradijs et 
al. 1995; Ghosh, Angelini, & White 1997). For all accretion 
models, there is an equilibrium pulsar spin period set by the 
mass accretion rate and the pulsar magnetic field strength. For 
the observed pulse periods and luminosities, the derived surface 
magnetic fields are of order 10 11 G. 

On the other hand, magnetar models have also been pro- 
posed, in which the AXPs are ultramagnetized (B~10 I4 -10 15 
G) neutron stars whose spin down is due to magnetic dipole 
radiation. In these models, the X-ray luminosity (which far ex- 
ceeds the spin-down energy) arises from magnetic field interac- 
tions in the neutron star crust (see, e.g., Thompson & Duncan 
1996). The spectral and timing properties of AXPs have been 
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used to constrain models for the atmospheres of ultramagne- 
tized neutron stars (Ozel, Psaltis, & Kaspi 2001). Other recent 
work on magnetar atmospheres has shown that proton cyclotron 
absorption features should be apparent in the soft X-ray spec- 
trum of magnetars (Ho & Lai 2001; Zane et al. 2001). These 
lines are predicted to be broad (AE/E w 1; Ho & Lai 2001). 
High-resolution soft X-ray spectroscopy should be able to test 
the magnetar model for AXPs. 

The brightest known AXP is 4U 0142+61 (/ = 129.°4, b = 
-0.°4), with a spin period of 8.7 s (see, Gavriil & Kaspi 2002, 
and references therein). Although long ago identified as an 
unusually soft X-ray source (Markert 1976; White & Mar- 
shall 1984; White et al. 1987), confusion with the nearby 
accretion-powered binary pulsar RX J0146. 9+6121 contami- 
nated the spectrum measured by early X-ray missions. More 
recent ASCA and BeppoSAX spectral studies measured a two 
component spectrum consisting of a 0.4 keV blackbody and a 
T=3.7 photon power-law (White et al. 1996; Israel et al. 1999; 
Paul et al. 2000). Hulleman, van Kerkwijk, & Kulkarni (2000) 
identified an optical counterpart based on the Einstein X-ray po- 
sition, and 8.7 s optical pulsations have recently been reported 
(Kern & Martin 2001). 

In this Letter, we discuss our Chandra X-ray Observatory 
High Energy Transmission Grating Spectrometer observation 
of 4U 0142+61. This is the first grating observation of an AXP 
obtained with Chandra. We show that the X-ray spectrum is 
consistent with previous lower resolution results and report on 
a search for line features in the spectrum. We place stringent 
limits on the strength of absorption lines and consequently the 
allowed magnetic field strengths assuming a magnetar model. 

2. OBSERVATION AND DATA REDUCTION 

We observed 4U 0142+61 with Chandra on 2001 May 23 for 
25 ks using the High Energy Transmission Grating Spectrome- 
ter (HETGS) and the spectroscopy array of the Advanced CCD 
Imaging Spectrometer (ACIS-S). For a more detailed descrip- 
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tion of the instruments, we defer to available Chandra X-ray 
Center (CXC) documents 3 . The HETGS employs two sets of 
transmission gratings: the Medium Energy Gratings (MEGs) 
with a range of 2.5-31 A (0.4-5.0 keV) and the High Energy 
Gratings (HEGs) with a range of 1.2-15 A (0.8-10.0 keV). The 
HETGS spectra were imaged by ACIS-S, an array of 6 CCD 
detectors. The HETGS/ACIS-S combination provides an undis- 
persed (zeroth order) image and dispersed spectra from the grat- 
ings. The various orders overlap and are sorted using the intrin- 
sic energy resolution of the ACIS CCDs. The first-order HEG 
(MEG) spectrum has a resolution of AA = 0.012 A (0.023 A). 
The ACIS CCDs are normally read out every 3.2 s. In order to 
sample the 8.7-s pulse period of 4U 0142+61 more effectively, 
we designed the observation to read out only a 512-row subar- 
ray on the ACIS-S detectors, resulting in a frame time of only 
1.8 s. 

The "level 1" event file was processed using the CIAO v2.1 
data analysis package 4 . During our analysis, it was found that 
the pipeline tool acis_detect_afterglow rejects 3-5% 
of source photons in grating spectra 5 . Afterglow is the residual 
charge left from a cosmic -ray event which is released over sev- 
eral frames and can cause a line-like feature in a grating spec- 
trum; the tool attempts to identify afterglow by flagging events 
that occur at the same chip coordinates in consecutive frames. 
However, for bright sources like 4U 0142+61, the tool also re- 
jects a small fraction of the valid source events. Although only 
a small fraction of the total, the rejection of source photons by 
this tool is systematic and non-uniform. Since order-sorting 
of grating spectra provides efficient rejection of background 
events, the afterglow detection tool is not necessary. Therefore, 
we reextracted the event file, retaining those events previously 
tagged by the acis_detect_afterglow tool in order to 
improve statistics. No spectral features were found that might 
be attributable to afterglow events. 

The zeroth-order image of 4U 0142+61 was affected by 
photon pileup (see, e.g., Davis 2001b) and was not used in 
our spectral analysis. However, we used the zeroth-order 
image to measure the source position, with the CIAO tool 
celldetect. Our best-fit position was R.A.=01 h 46 m 22?44 
and Dec. =6 1 45' 03 "3 (equinox J2000.0), with an approximate 
error radius of 0."5. The optical counterpart of Hulleman et al. 
(2000) lies 0."46 from our X-ray position. 

The combined MEG and HEG first order count rate was 
only 4.2 cts s" 1 . We therefore did not expect the dispersed 
spectrum to be affected by pileup, although we made sev- 
eral checks nevertheless. First, the dispersion distance, d, of 
each event was compared to the energy determined by the in- 
trinsic CCD resolution, £acis- The grating equation predicts 
£acis = hc/X = hcmR/(Pd), where R is the Rowland distance, 
m is the grating order, and P is the grating period. Sources 
with pileup show a deviation from this relationship in lower 
orders, as low order photons combine to produce higher en- 
ergy events which are then attributed to higher orders. There 
were no signs of pileup from the observed dispersion-energy 
relationship. Also, separate model fits of the first and second 
order spectra were found to be consistent both with past (low- 
resolution) observations and with each other, with no signatures 
of pileup in the residuals. We conclude that pileup is negligible 
in our dispersed spectra. 

3 http://asc.harvard.edu/udocs/docs/docs.html 

4 http://asc.harvard.edu/ciao/ 

5 http://asc.harvard.edu/ciao/threads/acisdetectafterglow/ 



We used the standard CIAO tools to create detector response 
files (ARFs; see Davis 2001a) for the MEG and HEG +1 and -1 
order spectra. These were combined when the +/- order spec- 
tra were added for the HEG and MEG separately. The spectra 
were binned in two different ways. For continuum analysis, we 
binned the data at 0.08 A with a minimum of 50 counts per 
bin. To look for high-resolution spectral features, the data were 
binned at 0.015 A for the HEG and 0.03 A for the MEG. We 
also created background files for the HEG and MEG spectra us- 
ing the standard CIAO tools. Finally, we used the XSPEC vl 1 
data analysis package (Arnaud 1996) to fit continuum models 
to the background-subtracted spectra. 
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FIG. 1. — (upper panel) Unfolded MEG and HEG energy spectra of 
4U 0142+61 fit with an absorbed power-law + blackbody model (solid line). 
The spectra have been coarsely binned for clarity. The contributions of the 
power-law (dashed line) and blackbody (dotted line) components are also 
shown. The normalization difference between the MEG and HEG is consis- 
tent with the uncertainty in the absolute flux calibration, (lower panel) Frac- 
tional residuals ([data-model]/model) from the absorbed power-law + black- 
body model fit. The feature at fs2 keV, marked with an arrow, is instrumental. 

3. SEARCH FOR ABSORPTION FEATURES 

We fit a variety of one and two component models to the 
coarsely binned HEG and MEG spectra: power-law, cutoff 
power-law, thermal bremsstrahlung, power-law + blackbody, 
power-law + disk blackbody, power-law + bremsstrahlung, 
bremsstrahlung + blackbody, and cutoff power-law + black- 
body. All models included absorption by neutral interstellar 
gas with cosmic abundances. The Chandra best fit parameters 
of the most acceptable models are given in Table 1 . 

The cutoff power-law + blackbody best fit parameters were 
inconsistent and gave a cutoff energy far above the instrument 
bandpass, suggesting that a standard power-law model was suf- 
ficient. The one component models had larger reduced \ 2 val- 
ues than the two component models and showed residuals that 
suggested the need for a second component. Using an F-test 
(see, e.g., Bevington & Robinson 1992), we determined that the 
addition of a blackbody to the power-law and bremsstrahlung 
models is highly significant. Therefore, we rule out the sin- 
gle component models: power-law, cutoff power-law, and ther- 
mal bremsstrahlung. The two component bremsstrahlung mod- 
els all gave large emission measures ss 10 58 c^ pc cm" 3 which 
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Table 1 



Continuum Spectral Fits' 





Nh 


T or kT 




Ecut 


kTbb 


Rbb 




Model b 


(10 22 crrT 2 ) 


( or keV) 


N c 


(keV) 


(keV) 


(km) 


X 2 /dof 


PL 


1.43±0.03 


4.02±0.05 


4.7+0.2 








515/259 


PL + BB 


0.88±0.13 


3.3±0.4 


1.0±0.5 




0.418±0.013 


2.0±0.2 


230/257 


BREMS 


0.92+0.02 


1.11 ±0.02 


1.13±0.05 








250/259 


BREMS + BB 


0.69±0.07 


1.6±0.4 


0.3±0.2 




0.391±0.014 


2.2±0.3 


227/257 


CUTOFFPL 


0.83±0.06 


0.9±0.3 


3.70±0.19 


0.93±0.11 






244/258 



J A11 errors are quoted at the 90%-confidence level 

b PL: power-law, BB: blackbody, BREMS: bremsstrahlung, CUTOFFPL: cutoff power-law 

L Power-law normalization at 1 keV in units of 1CT 1 photons keV -1 cnT 2 s -1 , bremsstrahlung normalization=emission measure 
in units of 10 58 cm 4 . 



would require an emission radius orders of magnitude larger 
than a neutron star for an optically thin plasma (White et al. 
1996), and were rejected on this basis. Traditionally, the source 
has been fit with an absorbed power-law + blackbody model 
and our analysis favors this model as well. The results of the 
Chandra fit, summarized in Table 1 and shown in Figure 1, are 
consistent with previous observations of 4U 0142+61. When 
the same model is applied to the high-resolution binned data 
set, the same parameter values are found to within the uncer- 
tainties. 



MEG 



0.008 

0.006 

< 0.004 

m 0.002 
0.000 

I 4 

c 0.008 
o 

° 0.006 

x 0.004 

z; 

c 0.002 
0.000 



l+ttltn 



' 2 



HEG 



10 



12 



Wavelength (A) 



FIG. 2.— MEG and HEG spectra of 4U 0142+61 at 6 X and 12 X instrument 
resolution, respectively. The best fit absorbed power-law + blackbody model is 
shown as a dashed line. The misfit of the continuum between 5.5-7.0 A is not 
consistent in the MEG and HEG spectra, and is attributed to instrumental fea- 
tures around known absorption edges. The spectra show no deep features that 
might result from proton cyclotron resonance absorption in a magnetic field in 
the range of (1.9-9.8)xl0 14 G. 



different depth in the MEG and HEG spectra (see Figure 2). 
For these reasons, we conclude that the feature is instrumental. 

To search for broad features, we used the high-resolution 
HETGS spectra and fit Gaussian line models to the fractional 
residuals ([data-model]/model) of the continuum fit. The cen- 
tral energies and widths of the Gaussian components were 
fixed, while the normalizations were fitted. To look for fea- 
tures that had been predicted in the magnetar models (Ho & 
Lai 2001; Zane et al. 2001), the sigma of the Gaussian was 
chosen to vary with energy {a = 0.1 x E). The Gaussian model 
was fit to the data centering at every wavelength point. We were 
able to determine the best fit amplitude and standard deviation 
of this result as well as the significance of each feature. The 4er 
upper limits on the equivalent widths (EWs) of any features are 
shown in Figure 3. There were no features with a significance 
greater than 4er, which we believe is a reasonable lower limit. 
In Figure 3, we have also included lines which mark the equiv- 
alent width limits for EW=0.1£ and EW=0.5£ for comparison. 
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With this continuum model, we then looked for absorption 
and emission features in the high-resolution spectrum. We 
searched for both narrow and broad line-like features. An ini- 
tial inspection of the residuals of the power-law + blackbody 
model fit seemed to show a feature at «6 A (2 keV; see Figure 
1). This feature is however highly suspect, since it is located 
between two known instrumental edges (Si-K at 6.75 A and 
Ir-M at 5.98 A). It does not appear in the MEG +1 spectrum 
(which is on a back-illuminated chip) but only in the MEG -1 
and HEG +1/-1 spectra (on front-illuminated chips). Further- 
more, the dip occurs at a slightly different energy and with a 



FIG. 3. — Upper bounds on the equivalent widths of lines with intrinsic 
width a = 0. IE (where E is the line energy). Magnetar atmosphere models pre- 
dict that the equivalent width of a proton cyclotron absorption line is of order 
0.70-0.75£ (Zane et al. 2001), although more recent work predicts equivalent 
widths at least an order of magnitude lower (Ho & Lai 2002). For comparison, 
we have included lines which define EW=0.5£ and EW=().1£. We conclude 
that the spectra show no strong absorption feature attributable to the proton 
cyclotron resonance in a magnetic field in the range of (1.9-9.8)x 10 14 G. 

Similarly, no narrow line features were found in the high- 
resolution spectra. We can place a 3a upper limit of 10 eV 
on any emission line between 4. 1-17.7 A (0.7-3.0 keV) whose 
width (FWHM) is comparable to the instrument resolution (see 
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Figure 4). The computation uses the model fit to the time- 
averaged HEG and MEG data. We assume that candidate fea- 
tures are only 2 bins wide. One may derive limits on broader 
features until the scale of the feature becomes comparable to 
that of the instrument calibration uncertainties. Limits on ab- 
sorption lines are identical when there are >25 counts per bin 
but are systematically larger for <25 counts per bin. 
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FIG. 4. — Upper bounds on the equivalent widths of narrow lines with 
FWHMs comparable to the instrument resolution, 0.023 A (0.012 A), for MEG 
(HEG). We assumed that candidate features are only 2 bins wide. One may de- 
rive limits on broad features using these curves until the scale of the feature 
becomes comparable to that of the instrument calibration uncertainties. Limits 
on absorption lines are identical when there are >25 counts per bin but are 
systematically larger for <25 counts per bin. 

4. DISCUSSION 

The Chandra spectrum of 4U 0142+61 is well fit by an ab- 
sorbed power-law + blackbody, which is consistent with pre- 
vious observations of the source (White et al. 1996; Israel et 
al. 1999; Paul et al. 2000). The small radius of the black- 
body model might be explained by atmospheric effects which 
give an overestimation of the temperature and thus an underes- 
timation of the radius when fit by a standard blackbody (see, 
e.g., Perna et al. 2001, and references therein). The best fit hy- 
drogen column density of (0.88±0.13)x 10 21 cm" 2 is consistent 
with estimates derived from radio and optical studies (Dickey 
& Lockman 1990; Hulleman et al. 2000). Distance estimates 
from extinction measurements for 4U 0142+61 range from 1-5 
kpc due to the line of sight passing through the edge of a lo- 
cal (< 1 kpc) molecular cloud. When including effects of the 
molecular cloud, a lower limit of 1 kpc is determined, whereas 
when neglecting the molecular cloud a lower limit of 2.7 kpc is 
obtained (Ozel et al. 2001). Normalizing the distance to 1 kpc, 



we find a luminosity (0.5-10 keV) of 5.3 xlO 34 dl pc ergs s 1 
which is consistent with past measurements from ASCA and 
BeppoSAX. 

No significant discrete features were found in the high- 
resolution spectrum. We can place a limit on the equivalent 
width of any narrow feature of «10 eV in the range 4.1-17.7 A 
(0.7-3.0 keV) and «50 eV for broad features in the range 2.5- 
13 A (0.95-5.0 keV). The limits on the equivalent width of 
any absorption features place strong constraints on the neu- 
tron star atmosphere models. Atmosphere models of ultram- 
agnetized neutron stars predict a broad proton cyclotron ab- 
sorption feature at an energy £ B = 0.63y i; (B/10 14 G) keV, where 
y g = (1 -2GM / c 2 R) l l 2 is the gravitational redshift factor and B 
is the surface magnetic field strength (Ho & Lai 2001; Zane et 
al. 2001). If we assume that magnetic dipole radiation drives 
the spindown, the magnetic field strength found using P and P 
(B 2 oc PP) is 1 .3 x 10 14 G (Gavriil & Kaspi 2002). A dipole field 
produces a cyclotron absorption feature consistent with a sur- 
face magnetic field strength 20-30% lower (Zane et al. 2001). 
For 4U 0142+61, the dipole field of 1.3 x 10 14 G would give a 
cyclotron absorption feature at «0.5 keV, which is well outside 
of the range accessible to our observation. 

Initial theoretical work predicted a proton cyclotron absorp- 
tion feature of equivalent width 0.10-0.1 5E B (Ho & Lai 2001; 
Zane et al. 2001), which is much greater than the equivalent 
width limits in the 0.95-5.0 keV range. If we take y g ~ 0.8, the 
above energy range corresponds to magnetic field strengths of 
(1.9-9.8)x 10 14 G. More recent work has shown that vacuum 
polarization effects strongly suppress the cyclotron absorption 
feature, giving equivalent widths an order of magnitude lower 
(Ho & Lai 2002). It should also be noted that the theoretical cal- 
culations of the equivalent width of the cyclotron feature were 
done for a local patch of the neutron star surface. A phase av- 
eraged spectrum, like the one we are using, would include con- 
tributions from various magnetic field strengths, directions and 
effective temperatures, which would further suppress the cy- 
clotron feature. The combination of these effects may reduce 
the predicted equivalent widths to below our measured values 
for much of the spectral range we studied. 
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tracts SAO SV1-61010 and NAS8-38249 as well as NASA 
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